Coupled layered superconductor as a system of 2D Coulomb particles of two kinds. 
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It is shown that the Josephson subsystem of the Lawrence-Doniach model of layered supercon- 
ductors in the London approximation can be presented as a system with variable number of classical 
Coulomb particles. This allows us to consider the vortex system of a coupled layered superconductor 
as the system of these particles and 2D-vortices interacting with each other. The grand partition 
function of the system was written and transformed into the form of field one. Thermodynamical 
properties of the model obtained was studied. It is found that there is no a phase transition in the 
system. Instead of this the model demonstrates the crossover from a low temperature 3D behavior 
to high temperature 2D one which can look as a phase transition for experimental purposes. 

PACS numbers: 74.20.-z, 74.25 Bt, 74.72.-h,64.60 Ak 



I. INTRODUCTION 

Behavior of the vortex system of coupled layered su- 
perconductors (SC) is one of the most interesting and 
difficult problem in the thermodynamics of quasi two di- 
mensional (2D) systems. It demonstrates dimensional 
crossover and other issues of dimensionality which are of 
perpetual interest in statistical mechanics. 

The Berezinskii-Kosterlitz-Thouless (BKT) phase 
transition^ in the 2D-vortex system was proposed to 
explain results of the experiment on measuring of the 
current-voltage characteristics of BISCO monocrystal^. 
Such a type of transition takes place in 2D systems in 
which topological defects can exist. They can be consid- 
ered as 2D Coulomb particles whose interaction energy 
logarithmically depends on a distance between them. 
Layered SC's without Josephson coupling between lay- 
ers are 2D systems. 2D-vortices are topological defects 
in these systems. 2D vortex-antivortex dipoles of a finite 
size can arise as thermal excitations at any temperatures 
and they dissociate into a gas of free vortices above BKT 
transition temperature T > Tkt- Properties of these 
systems are investigated very well^. 

But in real layered SC's layers are coupled due to the 
Josephson effect. They can not be considered as 2D sys- 
tems. Simplest thermal excitations in coupled layered 
SC's arise in the form of two 2D-vortices of opposite ori- 
entation localized in the same layer which magnetic flux 
is closed by two Josephson vortices or fluxons. The inter- 
action energy of 2D-vortices in this case logarithmically 
depends on a distance for a small dipole size and asymp- 
totically tends to linear dependence for a large size^. In- 
vestigation of such systems is much more difficult and 
less developed. 

One of the way which was used to analyze the problem 
consists in studying the thermodynamic of the system 
of classical particles with the pair interaction potential 
containing both logarithmic and linear terms (see for ex- 
ample Refs. 6.7,8]). This approach takes into account 
the linear tension of Josephson vortices as a portion of 
a two-particle interaction energy only while to describe 
the problem correctly it is necessary to take into consid- 



eration a Josephson vortex subsystem as an independent 
participant of thermodynamic processes which has own 
configuration energy, entropy, self-energy and interacts 
with the 2D vortex subsystem. 

The other approach to the problem is to investigate 
the anisotropic 3D XY-model. It was shown that excita- 
tions of the same type as in coupled layered SC can exist 
in this mode ffi 10 ' 11 . The anisotropic 3D XY-model was 
studied numerically using Monte Carlo simulation s 1 ^ 12 
and analyticall y 13 i 14 . Results of this investigations show 
that the model undergoes the phase transition of the 
BKT type at temperature T c higher then Tkt of a pla- 
nar model and coupling of adjacent layers vanishes at the 
same temperature. 

The attempt to solve the same problem, which our 
work is devoted to, was made in Ref. [HI]. The authors 
base their approach on the quantum ID sine-Gordon 
model and use the connection of this model with the 2D 
Coulomb gas system. The partition function of the model 
is represented in terms of a functional integral over two 
non-commuting variables. The model behavior is studied 
by means of a perturbative renormalization group (RG) 
approach. The RG recursion relations (RR) derived in 
the work have not a fixed point. However, authors an- 
alyzing RR features come to conclusion that the system 
undergoes phase transition which temperature Td is con- 
trolled by 2D-vortex core energy. 

The other approach to the problem discussed was pro- 
posed by Pierson and Valls 1 ^. The authors considered 
the model which they called a model of XY layers with 
Lawrence-Doniach (L-D) type interlayer couplings. It 
consists of 2D sine-Gordon models in each layers and 
additional cosine terms with difference of phases of ad- 
jacent layers, which describe the Josephson coupling. To 
study the model behavior authors used RG approach pre- 
viously expanding the additional cosine term and keeping 
the quadratic summand only. The RG RR's have not a 
fixed point at a small value of the coupling constant but 
authors give arguments for existing of that at a larger 
value beyond the range of validity of RR's. 

From our point of view such an approach to analyze 
the model is inadequate to the problem discussed. As 
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Nandori et ali^i showed in the case of the quadratic in- 
tcrlayer coupling the model is equivalent to the system 
of classical particles with long-range pair interaction. So, 
such an approach excludes from consideration the contri- 
bution of the entropy of the Josephson subsystem and the 
linear tension in the 2D-vortex interaction energy. 

In this work we show that the model considered by 
Pierson and Vails 16 follows from the L-D models in the 
London approximation and that its properties is deter- 
mined by behavior of order parameter singular points of 
two kinds. They can be considered as classical Coulomb 
particles which can be associated with 2D-vortices and 
the Josephson subsystem. 

The paper is organized as follow. In Sec. II we formu- 
late the L-D model in the London approximation which 
describes the system of Josephson vortices and transform 
it into the form available for further generalization. In 
Sec. Ill we introduce 2D-vortices in the L-D model and 
show the close analogy between the 2D-vortex and the 
Josephson systems. In Sec. IV the model of the vor- 
tex system of a coupled layered SC is constructed as a 
generalization of models considered in Sec. II and III. 
The behavior of the model is analyzed in Sec. V and 
VI by means of the RG approach and in the mean field 
(MF) approximation. Concluding remarks are presented 



in Sec. VII. 



II. FLUXON SYSTEM IN THE 
LAWRENCE-DONIACH MODEL. 

A. Model transformation. 

We start from the L-D model- in the London approxi- 
mation in which fluctuations of the order parameter mod- 
ulus responsible for transition in a superconducting state 
are neglected. In this case the order parameter is nor- 
malized so that its modulus is unit and does not depend 
on coordinates. Such a model describes only the sys- 
tem of Josephson vortices in a layered superconductor 
because 2D-vortices are associated with singular points 
in which the order parameter turns into zero. In this sec- 
tion we transform the model into the form suitable for 
its generalization and discuss briefly its thermodynamic 
properties. 

The Hamiltonian is the functional of the phase of the 
order parameter and the vector-potential A of magnetic 
field 



H = 



16tt 3 A(T) 



drdz 



V6»„(r) - x"A(r, z) I - 



-^2-cos(ft„(r)) 



5(z- 



dvdz 



[V x A(r,z)f 
8^ 



(1) 



Here cf>o is the magnetic flux quantum, A(T) = 
2A ab (T)/s, Aj = s\ c /X abl X a b(T) and A C (T) is the Lon- 
don penetration depths parallel and perpendicular to lay- 
ers, s is the period of layered structure and fl n (r) = 

0„+i(r) - e n (r) - (2vr/0 o ) J^ +1)S dzA z (r, z) is the gauge 
invariant phase difference. The partition function of the 
model can be written in terms of a functional integral as 



Z = J D6DAexp{-pH[0,A]} , 



(2) 



where (3 = 1/T. The vector potential A in this expres- 
sion is a Gaussian variable an can be integrated out. The 
variable A is a gauge field and this must be taken into 
account on calculating. In the case considered it is con- 
venient to choose the gauge condition A z = 0. Resulting 
partition function is 



W„(r)K„,„,VMr) 



Z = J D0expJ-i£ J dr 

| n,n' 

+ 2y f J2j Jcos(0„ +1 (r)-0„(r)) 



(3) 



where J = 47t 2 AT/</>q, yf is the fluxon fugacity which 
is proportional to Josephson critical current density and 
the operator 



K i = ft i — 



71 dk cos (k(n — n')) 



2tt 



1 + ±k 2 



takes into account magnetic interaction between layers. 
In systems we are interested the condition s <C X a b is 
fulfilled and the magnetic interaction leads to small cor- 
rections of order s 2 /A 2 and does not change results qual- 
itatively. Therefore for the sake of simplicity we will 
omit all over the paper these corrections and will take 
K n n' — &nri ■ For convenience of further references we 
write the partition function of this reduced model be- 
cause it is the initial point of our constructions 

+ / §cos(0„ +1 (r)-0 n (r)) j (4) 

n ) 

The field partition function ^ can be transformed into 
the form of that of a gas of classical Coulomb particles. 
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In order to perform such a mapping we follow to the 
scenario used by Nandori et al^i. With this aim in view 
we carry out the change of the variable 9 — \fjtp and 



J 



present the exponential function with the cosine term in 
its argument as a product of two exponents 



exp 



dr 



VJ((p n +i -<p n ) \ = fjexp 



dr 



,%V7(tp„+l—tp n ) 



yf t _ e v ™ ™ f exp^ V/ , nr 



dr 



Then we expand them into Tailor series. In the result after obvious change of summation we obtain the expression 
required 



n 



£ 

M„± 



1 1 

M n+ \M n A 



M n+ +M n 



J Dipexp | -- £ / dr 



2tt 



n j„=l 



r 



(5) 



In this expression g Jrl are charges of the fluxon 

particles, M n+ (M n _) are numbers of particles with a pos- 
itive (negative) charge in nth layer, /x/ is the chemical 
potential of the particles and length a is replaced by the 
correlation length £. 

The particles under consideration, apparently, are not 
real ones. They look as dipoles which poles have the same 
in-plain positions but are located in adjacent layers and 
can be presented as a point current source and sink. In 
fact, such a particle is a piece of a current line connecting 
adjacent layers. 

Further transformation of the partition function is pos- 
sible if we notice that ip is the Gaussian variable and can 
be integrated out. Such a procedure leads to the parti- 
tion function of the system of classical particles with a 
pair interaction 



in which the conditions 



n 



E 



1 1 

M„+! M n J. 



exp ] ~^E E 




The potential v nn > corresponds to 2D Coulomb inter- 
action of above-mentioned fluxon particles 



v nn ' (rj 



hi- 



Jn,n'-\-l 



■ (7) 



There is another way to obtain the partition function 
©■ We can start from Eq. ^ with yf = and as- 
sume that the order parameter possesses singular points 



(V,VMr-r,J) 
[V,V0 n (r)] 



±2tt JS(r - 




n-\-l,n f 



(8) 



are satisfied. It is easy to show that the energy of in- 
teraction of singular points defined by these conditions 
is the same Coulomb potential ([7]) as that of particles 
discussed. 

Thus, the field model ([3]) of fluxon subsystems can be 
presented as the system with variable number of classi- 
cal particles with pair interaction. Such a point of view 
on the fluxon system will be very useful in Sec. IV to 
construct the model in which fluxon and 2D-vortex sub- 
systems are joined. 



B. Thermodynamic properties of the model. 

Peculiarities of behavior of the model will reveal them- 
selves in the properties of the joined model. To better 
understand them we discuss briefly the thermodynamic 
properties of the fluxon subsystem. The main interest 
will be directed to an order and temperature of the phase 
transition. Similar model was studied by Horovitai^. He 
showed that the second order phase transition of the 
BKT type takes place in the system. Here we discuss 
results of our examination of the model which are very 
close to Horovitz ones. 

We perform the momentum space RG study of Eq.(|4]). 
RR's are a the set of two equations in parameters of the 
initial Hamiltonian 



d Vf 
dr 
dJ 

~d7 



(2-J)Vf, 
-8ir 2 J 3 y 2 f . 



(9) 
(10) 
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The phase portrait of the set is plotted in Fig. 1(a). It 
is the conventional picture of flows in the vicinity of a 
fixed point of the saddle type which is defined by the 
conditions (J = 2, t// = 0). Existence of the fixed point 
means that a second order phase transition takes place in 
the system at the temperature defined by the expression 



T 



S 



(11) 



which follows from the condition J = 2. This tempera- 
ture is placed above than that of the BKT transition Tkt 
in the 2D-vortex system of decoupled layered SC and is 
much more close to the temperature T c of the transition 
of a sample into a superconducting state. Another fea- 
ture of the transition is the directions of the RG flows. 
They are shown by arrows in Fig. [TJa) and, as it will 
be shown in next section (see Fig. Ufa)), are opposite to 
those of 2D-vortex system. 

The fugacity yf turns to zero in the fixed point. As yf 
is proportional to Josephson critical current density this 
means that it vanishes in the transition point T = Tt too 
and superconducting layers becomes decoupled at high 
temperatures T > Tf. 




0.94 0.96 
T/T 



FIG. 1: (a) Phase trajectories of the fluxon system RG equa- 
tions in the vicinity of the fixed point (J = 2, yf — 0) of 
saddle type, (b) Temperature dependence of the fluxon par- 
ticle concentration for the systems with different anisotropy 
c/ = 50 and 20 not far from the phase transition 

Let us consider also the temperature dependence of 
the concentration of fluxon particles. We found it in 



the MF approximation. The free energy of the system 
per the unit area of one layer was obtained from Eq.® 
by means of the thermodynamic perturbation theory in 
the ring approximation 20 . Under the suppositions that 
the system is neutral, equilibrium numbers of particles 
doesn't depend on the layer index (M n+ = M n _ = M) 
and the screening length Sf 3> £ the expression for the 
dimensionless free energy takes the form 



/ 



S 



2m (rn(m) — 1) — Jmln(4Jrn) + mJef. 



Here m = tt^ 2 M/ S is the dimensionless particle concen- 
tration, S is the layer area, and 1/yAJrn = Sf/£ is the 
dimensionless Debye screening length. The last term is, 
in fact, the chemical potential one. It controls the system 
anisotropy and we consider e/ as a phenomenological pa- 
rameter. The free energy (|12[) reaches its minimal value 
at the particle concentration 



4J 

e e f-> 



(13) 



This solution is stable at J < 2 (T < Tf). In the point 
J = 2 it loses the stability and at J > 2 (T > Tf) 
equilibrium value of the concentration is m — 0. So, the 
MF approximation leads to the same temperature (fTTj) of 
the phase transition in correspondence with RG results. 

The temperature dependence of the equilibrium par- 
ticle concentration is plotted in Fig. 1(b) for systems 
with different anisotropy. To derive this picture the 
model temperature dependence of the London length 
X a t(T) ~ (1 — T/T c ) _1 / 2 was used and parameters of the 
model were chosen in such a way to fix Tkt/T c = 0.98. 
These options result in Tf/T c « 0.998. It is seen that the 
free particle concentration m in a more anisotropic sys- 
tem (e/ = 50) always is less then that in a less anisotropic 
(e/ = 20). We can see also some correlations between 
temperature dependencies of the free fluxon particle con- 
centration and critical density of the Josephson current: 
both of them reach their maximum values at T = 0, both 
monotonically decrease when the temperature increases, 
both are greater in less anisotropic systems, both turn 
into zero at T — Tf and vanish at T > Tf . 



III. 2D- VORTEX SYSTEM IN THE 
LAWRENCE-DONIACH MODEL. 

A. Model constructing and transformation 

To describe the 2D- vortex system in the L-D model we 
use the Eq.(g]) with yf — 0. 2D-vortices are topological 
excitations and can be introduced in the model as singu- 
lar points of the order parameter placed in positions Tj n 
of the nth layer and defined by the conditions 

[V,V0 n (r)] = 2ns jn z5(r-r jn ), (V,V0„(r)) = 0. (14) 
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The symbol Sj n = ±1 means the sign of a charge of the 
vortex placed in the position Ta of the nth plane and z is 
the unit vector in z direction. The energy of interaction 
of two vortices defined in such a way and placed in the 
points Ti n and Tj n of the same nth layer is 



U(Ti 



-Qin Q] n U ( r i 



Q in Q jn In 



'Jr. 

- r,- 



(15) 



where Qj n = Sj n /\/~J is the 2D- vortex charge. The inter- 
action potential (fT5")l is the 2D Coulomb one and obeys 



the 2D Poisson equation 



Au(r) = 27T(5(r). 



(16) 



One can write the partition function of the system of 
2D-vortices considering them as classical massless parti- 
cles with the interaction potential d 



E 

N„+ 



1 



1 



N n+ \ N n _ 



N n+ +N n - 



exp 



,Qj n u(ri n - Tj 



(17) 



r 



Here N n ± are the numbers of particles with the posi- 
tive and negative sign of charge and fi v is the 2D-vortex 
chemical potential. This partition function can be trans- 
formed into the forms of that of a system of particles into 
self-consistent field and that of a field system^!. With 



this aim in view it is need to multiply the exponential 
function in Eq. p"?]) by the unit, which is represented as 
the ratio of two identical Gaussian integrals over a field 
variable <p> 



J 



exp 



^EE 



n i n ]„ 



Z 



Dtp exp 




(18) 



The integral in the denominator is just a normalizing 
constant Zq. The next step is to change of the variable 
in the numerator 



3n 



which is chosen in such a way that after integrating by 



J 



parts over r and using Eq. (|16[) it leads to compensation 
of the sum of interaction potentials in the power of expo- 
nential function (fT7)l . As a result, the partition function 
(fT7|) can be rewritten in the form of that of a system of 
particles in self-consistent field 



n 



E 

N„± 



1 



1 



dr 3u 

e 



Dip exp 



E 



2tt 



:>,< 



In the Eq. (|19p we can take the summation over particle numbers N n + and N n ~ to obtain the field sine-Gordon model. 
The substitution of the variable in the functional integration (p = dj\f~J returns us to the initial variable 6, which is 
the phase of the order parameter. The partition function finally takes the form 



6 



Here y v = exp{/3/z„} is the 2D-vortex fugacity and 
nonessential constant Z$ is omitted. 

Thus, the 2D-vortex system of a layered superconduc- 
tor can be represented in the same three forms as the 
fluxon one. Below we discuss briefly thermodynamical 
properties of the 2D- vortex system to compare them with 
those of the fluxon one. 



B. Thermodynamic properties of the model. 

Thermodynamical properties of this model was studied 
very much and are well knowni^. 

The momentum space RG approach based on Eq. ([2TJ1) 
gives the RR's 



dy v 
dr 
dJ 



1 

2J 



4tt^ 



(20) 
(21) 



Existence of the fixed point (J — 1/4, y v = 0) of the 
saddle type means that the second order phase transition 
takes place in the system at the temperature 



Trt = 



167r 2 A(T KT ) ■ 



(22) 



This is the conventional BKT transition. It is easy to see 
comparing Eqs. dHJ) and (11]) that T KT < T f . The phase 
portrait of the Eqs. (J20j) and ([21]) is plotted in Fig. 2(a). 
Comparison of Figs. 1(a) and 2(a) shows that RG flows 
of the fluxon and the 2D-vortex systems are oppositely 
directed. 

To find temperature dependence of the 2D- vortex con- 
centration we used the MF approximation. The free en- 
ergy of the system as a function of the vortex number 
can be obtained from the partition function Eq. (I17|) by 
means of the thermodynamic perturbation theory in the 
ring approximation 20 . Under supposition that the sys- 
tem is neutral the dimcnsionless free energy per the area 
unit 7r£ 2 of one layer takes the form 

/ = 2 „( h „ -!)-£(,„! _l) +B i. (23) 

Here n — ir£ 2 N/S is the dimensionless vortex concentra- 
tion, y J /An = 5 v /t; is the dimensionless Debye screening 
length, e v /4J is the vortex core energy. 

The equilibrium vortex concentration obeys the equa- 
tion 



— e 

J 



(24) 



This solution is stable at T > T KT . At T = T KT it 
becomes zero and loses the stability at T < Tkt where 
zero solution becomes stable. Temperature dependence 
of the concentration is shown in Fig. (Hb). To obtain this 




0.996 



FIG. 2: (a) Phase trajectories of the 2D vortex system RG 
equations in vicinity of the fixed point ( J = 1/4, y v = 0). (b) 
Temperature dependence of the free 2D-vortex concentration. 



picture we used model temperature dependence X a b(T) ~ 
(1 — T/T c )~ 1//2 and chouse model parameters to specify 
Tkt/T c = 0.98. 

It is interesting to compare temperature dependencies 
of fluxon particles Fig. []Jb) and 2D-vortices Fig. [2Jb) 
concentrations. A state with a finite concentration of 2D- 
vortices is stable in the temperature range Tkt < T < T c 
which is located higher the phase transition temperature 
while that of fluxon particles is stable at lower tempera- 
tures < T < T f . 

So, we considered two basic models. We found their 
common features and revealed differences between them. 
By this we finish consideration of the previous results. 
We have all we need to construct the model of the vortex 
system of a coupled layered superconductor. This is the 
subject of next section. 



IV. MODEL OF INTERACTING 2D- VORTEX 
AND FLUXON SYSTEM. 

Now as a result of previous consideration it is easy to 
see close analogy of the fluxon and 2D- vortex subsystems 
of a layered SC. Each of them can be presented either as 
a held model of a sine-Gordon type or as a system of 
classical particles with a pair Coulomb interaction. Dif- 
ference between them consists in the structure and the 
value of "charge" of the particles only. But our aim is 
to develop the model describing a coupled layered SC 
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which contains both these subsystems. Therefore we will 
consider a coupled layered SC as a system consisting of 
two kinds of 2D Coulomb particles, namely, fluxon par- 
ticles and 2D-vortices, interacting with each other. Such 
a view on layered SC's allow us to obtain the partition 
function of their vortex system to study thermodynamic 
properties. 

There are two ways to solve this problem. One can 
act in the complete analogy with constructing of the 2D- 
vortex subsystem in Sec. Ill taking into account singular 
points of the order parameter defined by both conditions 



(J5J) and ([M]) to obtain 2D-vortex, fluxon particle and 2D- 
vortex - fluxon particle pair interaction potentials and 
write a partition function. 

The other way is to construct the partition function 
of the system of particles of both kinds in self-consistent 
field. Really, existence of the interaction between the 
particles means that all of them produce the same field 
and interact with it in accordance with their structures 
and charges. In the analogy with Eqs. (O and (fT9|) we 
can write 



E 



1 1 

aTTT 7V i 

iv n+- iv n— • 



N n+ +N n - 



E 



Mr, 



1 1 

M n+ \M n _l 



M n+ +M„ 



J { n J ' n j n 



'( r i»)+'EE? a » (fn+i(r an ) ~ Vn(r a J) \ • (25) 



Next step is to carry out the sums over all numbers 
of particles. These sums are easy to calculate because 
the result of integration over particle coordinates r Jre and 
v an do not depend on particle indices and sums are just 



Taylor series of exponential functions. As a result, the 
partition function after change of a variable ip = 6/yJ 
takes the form 



D9 exp ■ 



-E 

2J ^ 



dr 



2tt 



2 ^e/§ 



cos ( -9„{t) 



dr 

— cos(0„+i(r) - 0„(r)) 



(26) 



The expression obtained is the model we developed. It 
describes the fluxon subsystem, the 2D- vortex subsystem 
as well as interaction between them and is, in fact, the re- 
duced L-D model in which order parameter fluctuations is 
neglected. This model includes the initial models, which 
was the base to develop the model, as limiting cases. In 
the case y v — the model (|26|) reverts back to that of 
the fluxon system ([4]) and the condition yj — convert 
it into Eq. ([2H)l which is the partition function of the 
2D-vortex subsystem. 

The model considered coincides actually with that pro- 
posed by Pierson et alJ£. To investigate the model au- 
thors expanded the second cosine term and performed 
a RG analysis of the reduced model. Nandori et al.— 
showed that the reduced model is equivalent to a gas 
of topological defects with long-range interaction poten- 
tials. But as it was shown in this section the model [26] is 
equivalent to the gas of classical Coulomb particles of two 
kinds. In contrast to worki^ we study the model taking 
this into account. 



V. RENORMALIZATION GROPE ANALYSIS. 

We derived the perturbative RG equations in the pa- 
rameters of the Hamiltonian of the model (|2l)|) by means 
of the momentum space approach. Recursion relations 
we obtained 



dJ 
~d7 

dy v 
dr 
dyf 
dr 



4tt^ 
2 



y 2 v -8ir 2 J 3 y}, 



1 

2J 



yv, 



(27) 
(28) 
(29) 



takes into account renormalization of the parameters of 
the initial Hamiltonian only. It is easy to see that the 
set of equations inherits the structures of sets ([ MTU]) and 
(I20H21[) as well as the model (|26p inherits the structures 
of models ^ and (f2T)]) . Equations in fugacities y v and 
yf are just identical and the right-hand side (r.h.s.) of 
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Eq. |57]| is the sum of those of Eqs. © and ([20]). Such 
a structure of the equations results in that some details 
of the model behavior are similar to that of the fluxon or 
2D-vortex models but another ones are quite different. 

Behavior of the set considered is much more compli- 
cated. Main peculiar features of the set are that the 
number of independent dynamical variables rises from 2 
to 3 and the r.h.s. of Eq. ([2~T|) is not sign-definite. In the 
result the recursion relations ([2"T)) - ([2lI|) , contrary to that 
of the fluxon system ©-(Hi]) and the 2D-vortex one ([2"0]) - 
(f2Tj) . have not a fixed point and thus the system does not 
undergo a second order phase transition. 

Typical trajectories of the RG recursion relations are 
shown in Fig. [3]as projections on planes y v — J and y/ — J- 
The initial conditions are chosen in such a way to get 
the curves which projections behavior, in a large extent, 
are similar to those of corresponding trajectories of the 
fluxon and the 2D-vortex models. It is easy to see that 
the projection of trajectory marked by 1 in Fig. [3Ja) at 
small values of the scale variable r behaves analogous to 
the trajectory 1 in Fig. Ela). But at some finite value of r 
second term in r.h.s. of Eq. (|27p becomes dominated and 
the derivative of J changes the sign from plus to mines. 
It is a turning point of the trajectory. After this point 
the analogy mentioned disappears but instead analogy 
between the projection 1 in Fig. [3Jb) and the trajectory 
1 in Fig. Ufa) appears. Similar analogies can be find 
between the projections of trajectory 2 in Fig. [3] and 
the trajectories 2 in Figs. HJa) anddja). There are else 
two analogies between phase portraits of these systems. 
The projection 3 in Fig. [3^a) is similar to the trajectory 
3 in the phase portrait of the RG equations of the 2D- 
vortex system, and the projection 4 in Fig. [31(b) behaves 
analogously to the trajectory 3 in the phase portrait of 
the RG equations of the fluxon system. 

Such an analogy between phase portraits of the model 
considered and the models of independent 2D- vortex and 
fluxon systems allows us to believe that there are some 
similarity between behavior of these systems. But this 
analogy is observed in the space of parameters of Hamil- 
tonian and it does not allow us to draw a valid conclu- 
sion about behavior of the vortex system in the space of 
thermodynamical variables. Let us consider the phase 
portrait shown in Fig. [3] from another point of view. 

The RG flow can be characterized, in addition to forms 
of trajectories and direction of the system motion along 
them, also by a rate of the system motion. If the recur- 
sion relations have a fixed point the rate of motion along 
the trajectory is minimal in the position which is nearest 
to this point and asymptotically tends to zero when a 
trajectory approaches to it. Such a situation takes place 
in the RG flows of the fluxon Fig. HJa) and the 2D- vortex 
Fig. [2a) systems. 

In the case of the model proposed the situation is dif- 
ferent. The RG flow Fig. [3] does not have a fixed point. 
Instead, it has a slow point in which a rate of motion 
along a trajectory is lowest. This point can mimic a 
fixed point for experimental purposes. The properties of 
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FIG. 3: Projections of phase trajectories of the model RG 
equations. 

models which demonstrate similar behavior were investi- 
gated by Zumbach22. He termed such a situation as an 
" almost second order phase transition" . 

From such a point of view we can suppose that the 
model asymptotically behaves almost the free fluxon and 
the 2D- vortex systems in regions of low (T < Tkt) 
and high (T > Tf) temperature and demonstrates the 
crossover from low temperature 3D behavior to high 
temperature 2D one in vicinity of temperature T s cor- 
responding to a slow point. The turning points of trajec- 
tories of different kinds (1-4 in Fig© draw together in 
close vicinity of a slow point. This allow us to evaluate 
the position of this point. Direct numerical calculation 
of trajectories passing in the vicinity of slow point shows 
that such a point corresponds to J « 0.7. The crossover 
temperature which can be obtained from this relation 
obeys the equation 



But a perturbative RG is not a very good approach to 
find a slow point because the position of this point can be 
dependent on a cut-off procedure. So, the expression ([50]) 
can be considered as a rough evaluation of the crossover 
temperature only. The fact in which we can be sure is 
that T KT <T S <Tf or the same 1/4 < J s < 2. 

In completion of the section we note that a RG ap- 
proach is the examination of the hypothesis about a scal- 
ing invariance of the model. In the case considered the 
RG does not have a fixed point and, thus the hypothesis 
is not prove to be true, in contrast to the independent 
2D-vortex and fluxon models. In this situation we have 
to use a different approach to understand reasons of the 
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scaling invariance breaking and to investigate the model 
behavior in more details. This is a subject of the next 
section. 

VI. MEAN FIELD ANALYSIS. 

In this section we will consider the model as a two 
subsystems of classical particles with the long-range 



Coulomb interactions ([7]) and (fTB"]) interacting with each 
other. To obtain the free energy of the model in a MF ap- 
proximation we use the method developed in Ref. [2a |. 
This method corresponds to the ring approximation in 
the case of the Coulomb gasSS,. The system free en- 
ergy obtained is a function of concentrations of both 
2D-vortices (n) and fluxon particles (to). Under the 
assumption that the system is neutral (N n+ (M n+ ) = 
N n _(M n J) = N(M)) the free energy takes the form 



/ = 2n(lnn - 1) + 2m(lnm + 2mJ ) ln ( J J 



j+*rnj)+~J 1 J 1 + 4mJ + n£ + mJe f . (31) 



The condition of a minimum of the free energy is two equations in equilibrium concentrations 

1 



2 ln n ln 

J 



2 In to - 2 Jin 



J 



AraJ 



AmJ 



2 J 



0. 



n_ m +4mJ 
J^L — V J + Jej 



(32) 



The solutions of the equations obtained numerically are 
shown in Fig[5]as temperature dependence of the equilib- 
rium concentrations of the 2D-vortices and fluxon parti- 
cles by solid lines for more anisotropic system (e/ = 50) 
and by dotted lines for less anisotropic (e/ = 20) . These 
concentrations always are above the corresponding val- 
ues of the independent subsystems, which are plotted by 
the dashed lines, and in contrast to those are finite at 
any temperatures in the whole region < T < T c . 

This result is easy to understand. There is BKT phase 
transition in the independent systems of 2D-vortices and 
fluxon particles. The mechanism realizing the transition 
is follow. A single 2D-vortex (fluxon particle) in an infi- 
nite sample has infinite energy. As a result of this such a 
particle can not appear in the system as a thermal fluc- 
tuation at low temperature T < Tkt (high one T > Tf). 
But in a many body system the competition between the 
configuration energy and entropy terms in the free en- 
ergy (see Eqs. and f23]|) at T > T KT (T < T f ) 
leads to instability of a zero concentration state and to 
appearance of a finite concentration of free 2D-vortices 
(fluxon particles). These particles differ from single ones 
mentioned above because they have a finite energy due 
to Debye screening. Such free vortices (fluxon particles) 
are, in fact, quasi-particles and can appear in the system 
as thermal fluctuations. 

The model discussed describes systems of 2D-vortices 
and fluxon particles interacting with each other. In the 
case of Coulomb particles the interaction leads to their 
mutual screening. Thus, a finite concentration of 2D- 
vortices (fluxon particles) at low temperature T < Txt 
(high one T > Tf) exists due to screening of them by 
fluxon particles (2D- vortices). This supposition is easy to 
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FIG. 4: Solutions of equilibrium equations are shown by solid 
lines for the high anisotropic system e/ = 50 and by dotted 
lines for the low anisotropic one e/ = 20. Dash and dash and 
dot lines show asymptotical solutions Eqs. (|33[1 and (|34[) 



verify. The equilibrium equations (f3"2"]) can be simplified 
in the vicinity of T = (T = T c ) due to condition n <C m 
(to <C n). Asymptotical solutions of the equation take 
the form 



4J 



a = exp <i — (ln4jTO — e v ) 



(33) 
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in the vicinity of T = and 




close to T = T c . 

The solutions obtained make a clear sense. The con- 
centration to in (|33p coincides with that in the inde- 
pendent fluxon system (fT3")) . This is sequent of a very 
small concentration n of 2D- vortices at low temperatures 
T <C Tkt, which can be neglected in comparison with to. 
The concentration of 2D-vortices n in the limit consid- 
ered is just the Boltzmann expression n = exp{— E v /T}, 
where E v is the vortex energy containing two terms. Sec- 
ond term is the vortex core energy g v /4:J. First one is 
the energy of a 2D-vortex which diverges logarithmically 
in the case of a single vortex but is limited by the screen- 
ing length Sf due to a finite concentration m of fluxon 
particles, AmJ = £ 2 ''/Sf. The Eq. (f34|) has an analo- 
gous structure: the 2D-vortex concentration n coincides 
with that in the independent system Eq. (|2"4")1 and the 
concentration m is defined by the self-energy of a fluxon 
particle screened by a finite concentration of 2D- vortices, 

Asymptotical temperature dependencies of the concen- 
trations of the particles of both kinds are in a good agree- 
ment with the numerical solution. The concentrations of 
2D-vortices and fluxon particles of independent systems 
are plotted in Fig. 0] by dashed lines. Asymptotical de- 
pendencies of concentrations of 2D- vortices at low tem- 
peratures and fluxon particles at high temperatures are 
shown by dash and dot lines. 

It is seen that asymptotical behavior of the model at 
T — > and T — > T c is very close to that of the indepen- 
dent 2D-vortex and fluxon systems. The concentrations 
natT< Tkt and m at T > Tf are exponentially small 
but finite. This means that there is no a phase transition 
in the model discussed but there is a crossover from 3D to 
2D type of behavior which takes place in the temperature 
interval Tkt < T < Tf . 

Thus results obtained in the framework of the MF ap- 
proximation agree quantitatively with that of the RG 
analysis. In terms of the MF approximation they are 
conditioned by mutual screening of particles of two kinds. 



The RG considers this as breaking of the scaling invari- 
ance of the model. Now we can see that the reason which 
leads to this is existence of two competitive lengthes, 
namely, the screening lengthes S v and 6 / . 

VII. CONCLUSION 

We constructed and studied the model of the vortex 
system of coupled layered superconductors which is based 
on the L-D one. The main idea is to consider the Joseph- 
son and 2D-vortex subsystems as the systems of singu- 
lar points defined by the conditions © and (|14[>. Such 
an approach becomes possible because both singularities 
interacts according to the 2D Coulomb law. So, they 
can be interpreted as classical massless Coulomb parti- 
cles which are characterized by their charges and have 
different structures. The model partition function can 
be represented either as a grand partition function of the 
system with a variable number of particles of two kinds 
(|2"5]1 or as a that of the field model (gUJ) . 

The model was examined by means of the perturba- 
tive RG approach and the MF approximation. Results 
obtained by both methods agree qualitatively with each 
other. Both approaches show that there is no phase tran- 
sition in the model in the whole temperature interval 
< T < T c were the model is defined. But the model 
behaves asymptotically at T — » and T — > T c as in- 
dependent 2D-vortex and Josephson systems. Crossover 
from low-temperature 3D behavior to high-temperature 
2D one takes place in the interval Tkt <T <Tf. 

The effect of anisotropy can be cleared by means of 
comparing of temperature dependencies of free particles 
concentrations of systems with different value of the pa- 
rameter ef in Fig. |4j One can see that more anisotropic 
system (e/ = 50) behaves closer to that of indepen- 
dent 2D-vortex and fluxon particle systems than less 
anisotropic (e/ — 20). If the layered system is anisotropic 
enough it can mimic the BKT phase transition for exper- 
imental purposes. 
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